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Abstract 

We have performed analysis of the solar neutrino data including results from the 
SNO salt phase as well as the combined analysis of the solar and the KamLAND 
results. The best fit values of neutrino parameters are Am^ = 7.1 • 10~^ eV^, tan^ 9 = 
0.40 with the boron flux fs = 1.04. New SNO results strongly disfavor maximal 
mixing and the h-LMA region (Am^ > 10~^ eV^) which is accepted now at the 3a 
level. We find the 3a upper bounds: Am^ < 1.7 • 10~^ eV^ and tan^ 6* < 0.64, 
and the lower bound Am^ > 4.8 • 10~^ eV^. Non-zero 13- mixing does not change 
these results significantly. The present data determine quantitatively the physical 
picture of the solar neutrino conversion. At high energies relevant for SNO and Super- 
Kamiokande the deviation of the effective survival probability from the non-oscillatory 
value is about 10 - 14%. The oscillation effect contribution to this difference is about 
10% and the Earth regeneration about 3 - 4 %. At low energies {E < 1 MeV) the 
matter corrections to vacuum oscillation effect are below 5%. The predictions for the 
forthcoming measurements are given which include the spectral distortion and CC/NC 
ratio at SNO, the Day-Night asymmetry, the KamLAND spectrum and rate. 



1 Introduction 



The SNO-II (salt phase) results [1] have further confirmed correctness of both the neutrino 
fluxes predicted by the Standard Solar Model (SSM) [2] and picture of the solar neutrino 
conversion based on the MSW effect [3, 4]. Together with the results from the SNO phase- 
I [5], Homestake [6], SAGE [7], GALLEX [8], GNO [9] and Super-Kamiokande [10, 11] 
experiments as well as from the reactor experiment KamLAND [12] the latest SNO results 
lead to better determination of the oscillation parameters. This improvement allows one to 
make two important qualitative statements [1]: 

• the h-LMA region, which corresponds to Am^ > 10~^ eV^, is strongly disfavored 
(being accepted at the 3a level); 

• there is a substantial deviation of 1-2 mixing from the maximal one, and the latter is 
rejected by more than 5a standard deviations. 

These statements have been confirmed by followed analyzes [13, 14, 15, 16, 17, 18]. 

There are two key points of the SNO-II publication: (i) new measurement of the total 
active Boron neutrino flux with enhanced neutral current sensitivity, and (ii) data analysis 
without assumption of the undistorted neutrino spectrum. They lead to larger (than pre- 
salt) value of the neutrino flux measured by the neutral current (NC) reaction and smaller 
flux measured by the charged current (CC) reaction [1]. As a result, the ratio of neutrino 
fluxes [1]: 

CC 

— = 0.306 ± 0m6{stat) ± 0.024{syst) (1) 

IN 

turns out to be smaller than the value from pre-salt measurements: 0.346 + 0.048/ — 0.046 
[5]. It is, however, larger than the value which corresponds to the NC flux determined 
previously without constraint of the undistorted spectrum: CC/NC = 0.274 ± 0.073 [5]. 

The ratio (1) is also smaller than the expected one. In the previous paper [19] we have 
predicted that the central value of the ratio (which corresponds to the best fit point) and 
the 3a error bars equal: CC/NC = 0.33 -|- 0.15/ — 0.07. It was pointed out that values 
CC/NC< 0.35 exclude the h-LMA region. Furthermore, precise measurements of CC/NC 
should strengthen the upper bounds on mixing and Am^ [19]. 

Confronting the prediction for CC/NC with the experimental result (1) one can under- 
stood consequences of (1) immediately. According to fig. 10 of rcf. [19] (where we show the 
contours of constant CC/NC ratio in the Am^ — tan^ plot) the ratio CC/NC decreases 
with Am^ and tan^ 9. Therefore, with the SNO salt results the best fit point and the allowed 
region shift toward the smaller values of Am^ and tan^ 9. Maximal mixing and the h-LMA 
region should be further disfavored. This is indeed confirmed by the detailed studies [1] . 

In this paper we present our analysis of the solar neutrino data including the SNO-II 
fluxes (sec. 2). We perform a combined analysis of the solar neutrino and KamLAND 
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results (sec. 3). Possible effect of 1-3 mixing is considered in sec. 4. We show that after the 
SNO-II measurements, physics of the solar neutrino conversion is essentially determined. In 
particular, relative contributions of the adiabatic conversion in matter and the oscillation 
effect for different energies can be quantified (sec. 5). In sec. 6 the expected spectrum 
distortion at SNO is studied. In sec. 7 we give predictions for the future experiments. The 
conclusions are presented in sec. 8. 

2 Solar Neutrinos with the SNO salt-phase results 

In the analysis we use the following data set: 

• 3 total rates: (i) the Ar-production rate, Qat, from Homestake [6], (ii) the Ge— production 
rate, Qce, from SAGE [7] and (iii) the combined Ge— production rate from GALLEX [8] 

and GNO [9]; 

• 44 data points from the zenith-spectra measured by Super-Kamiokande during 1496 
days of operation [10, 11]; 

• 34 day-night spectral points from SNO [5] ; 

• 3 fluxes from the SNO salt phase [1] measured by the CC-, NC and ES- reactions. We 
treat correlations of these fluxes following prescription in [20] . 

Altogether the solar neutrino experiments provide us with 84 data points. 

All the solar neutrino fluxes, but the boron neutrino flux, are taken according to SSM 
BP2000 [2]. The boron neutrino flux measured in the units of the Standard Solar Model 
flux [2], Jb = Fb/F§^^^, is treated as a free parameter (here F§^^ = 5.05 ■ 10^ cm"^ s~^). 
For the /lep— neutrino flux we take fixed value F^ep — 9.3 x 10^ cm"^ s"^ [2, 21] . 

We use the same procedure of analysis as in our previous publication [22, 19]. In analysis 
of the solar neutrino data as well as in the combined analysis of the solar and KamLAND 
results we have three fit parameters: Am^, tan^^ and fs- 

We perform the test defining 



where Xrate^ Xsk^ Xsno-i Xsno-ii contributions from the total rates, the Super- 

Kamiokande zenith spectra, the SNO day and night spectra and the SNO fluxes from the 
salt phase correspondingly. The number of degrees of freedom is 84 - 3 = 81. 

The minimum, ■x^g^^{min)/d.o.f. = 67/81 is achieved at 




Xrate + XSK + XSNO-I + XsNO-IIi 



(2) 



Am^ = 6.31 • 10-^eV^ tan^ 9 = 0.39, fs = 1.063. 



(3) 
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From the pre-salt analysis we had: {Am^,t£in'^e) = (6.15 • 10-^0.41) [22]. So, the SNO-II 
data only slightly shifted the best fit point to smaller values of Am^ and tan^ 9. The boron 
neutrino fiux was smaller: J'b = 1.05. 

We construct the contours of constant confidence level in the Am^ — tan^ 6 plot (fig. 1) 
using the following procedure. We perform minimization of with respect to Jb for each 
point of the oscillation plane, thus getting xLn(^^^' ^^^^ Then the contours are defined 
by the condition x^(Am^, tan^ ^) = x^(mm) + Ax^, where Ax^ = 2.3, 4.61, 5.99, 9.21 and 
11.83 are taken for la, 90%, 95% and 99% C.L. and 3a. 

In contrast to the best fit point, the influence of the SNO-II data on the size of allowed 
region, especially for high confidence levels is much stronger. The data reduce substantially 
the allowed region from the parts of large Am^ and tan^ 6. The borders of 3a region are 
shifted as (4.5 ^ 1.6) • 10"^ eV^ for Am^, and (0.84 ^ 0.64) for tan^ 9. 

Shift of the la contour is much weaker. Projecting the la region from fig. 1 we find the 
intervals: 

Am^ = (3.8 - 10) • 10-^eV^ tan^ 9 = 0.325 - 0.475. (4) 
The lower bounds for Am^ and tan^ 9 are practically unchanged. 



3 Solar neutrinos and KamLAND 

The KamLAND data have been analyzed using the Poisson statistics. Minimizing 



i=l,13 



- Nf' + NfHn 



m 



th 



(5) 



we find the best fit point 

Am^ = 7.24 • 10-^eV^ tan^ 9 = 0.52. (6) 

Confronting the oscillation parameters determined independently from the solar neutrino 
experiments (3) and KamLAND (6) one can check CPT and search for new physics "beyond 
LMA". Since there is no significant change in the best fit point and la region for solar 
neutrinos, status of the CPT check is practically unchanged in comparison with the pre-salt 
analysis. The data are well consistent with the CPT conservation. As before, there is an 
overlap of the la regions of oscillation parameters found from the solar data and KamLAND 
results. The b.f. point (3) is at the border of 95% C.L. region allowed by KamLAND [12]. 

Assuming CPT conservation we have performed a combined fit of the solar neutrino 
data and KamLAND spectral results. We calculate the global x^ 

^global — Xsun ~^ XkL " (''') 
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There are 84 (solar) + 13 (KamLAND) data points - 3 free parameters = 94 d.o.f.. The 
absolute minimum, x^iobaii™''^)/^-^-!- — 73.4/94 is at 

Am^ = 7.13 • 10-^eV^ tan2^ = 0.40, /s = 1.038. (8) 

The SNO salt results have shifted the best global fit point to smaller Am^ and tan^ 6 in 
comparison with the results of previous analysis [19]: Am^ = 7.30 • 10~^eV^, tan^ 9 = 0.41, 
fB = 1.05. 

We construct the contours of constant confidence level in the Am^ — tan^ 6 plot (fig. 
2) using the same procedure as in sec. 2. According to fig. 2 the intervals of parameters 
obtained by projection of the la allowed region equal 

Am^ = (6.4 - 8.4) • lO'W, tan=^ 9 = 0.33 - 0.48. (9) 



Notice that inclusion of the KamLAND result does not change the 3cr upper bounds on 
Am^ and tan^^ (compare with fig. 1). However, KamLAND further strengthen the lower 
bound on Am^: 

Am^ > 4.8 • 10"^ eV^ {3a). (10) 

As a result of the SNO-II measurements, the h-LMA region, Am^(12 - 16) • 10~^ eV^, 
is strongly disfavored being accepted with respect to the global minimum (8) at about 3a 
level only. 



Prom fig. 2 we get the following upper bounds on mixing: 

tan^ 9 < < 



0.48, la 

0.55, 2a . (11) 
0.64, 3a 



Maximal mixing is excluded at 5.2a. These bounds follow from the solar neutrino data 
being practically unaffected by KamLAND. 

To check stability of our results we have performed analysis taking the boron neutrino 
fiux as it is predicted by the SSM: fs — 1, with the corresponding theoretical errors [2]. 
Results of the analysis with two free parameters: Am^ and tan^ 9 are shown in fig. 3. In 
the best fit point: 

Am^ = 7.36 ■ 10-^eV^ tan^ 9 = 0.41 (12) 

we have x^/ d.o.f. = 74.0/82. The allowed region is shifted to larger tan^ in comparison 
with the free Jb fit. Indeed, in the latter case we had fs > 1, so to compensate the decrease 
of /b the survival probability, which is proportional to sin^ 9, should increase. Also the 
contours of constant confidence level shift to larger mixings. The 3a upper bound becomes 
weaker: tan^ 9 < 0.628. 

Notice that the h-LMA region is rejected now at the 3a level. 
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4 Effect of 1-3 mixing 



As it was found in the pre-salt phase analysis, the 1-3 mixing improves the relative goodness 
of the fit in the h-LMA region (see, e.g., [19]). Can non-zero 1-3 mixing substantially change 
the results of sect. 3? 

Both for the KamLAND and for solar neutrinos the oscillations driven by Amjg are 
averaged out and signals are determined by the survival probabilities 

Pee = (1 - sin^ ^13)2P2 + siu^ ^13 ^{1-2 siu^ ^lajPa, (13) 

where P2 = P2(Ami2, ^12) is the two neutrino vacuum oscillation probability in the case of 
KamLAND and it is the matter conversion probability in the case of solar neutrinos. 

The effect of non-zero 1- 3 mixing is reduced to the overall suppression of the survival 
probability. According to recent analysis of the atmospheric neutrino data [23] the allowed 
region of oscillation parameters has shifted to smaller Am^a, and consequently, the upper 
bound on 1-3 mixing from the CHOOZ experiment [24] becomes weaker: sin^^ia = 0.067 
(3(t) [25]. For this value of sin^ ^13, the suppression can reach 13% in (13). The influence of 
sin^ ^13 on the global fit can be traced in the following way. 

There are three sets of observables for which effects of sin^ ^13 are different. 

1) Total rates (fluxes) at high energies measured by SNO and SK. These rates depend 
essentially on the combination cos^^i3P2. In particular, 

CC ES 

— ^ cos" ^13P2, ^ ~ COS^ ^13^2(1 - r) + T. (14) 

The ratios of fluxes are unchanged if 

cos^ ^i3(P2(Ami2, tan^ ^12)) = const., (15) 

where (...) is the averaging over the relevant energy range. The product (15) is invariant 

with increase of sin^ ^13 (decrease of cos^ ^13) if the probability increases. In the region near 
the b.f. point (8) the latter requires increase of Am^2 or/and tan^^i2. Then the absolute 
values of fluxes can be reproduced by tunning fs (in the free Jb flt). 

With increase of sin^ ^13 the predicted spectrum becomes flatter which does not change 
the quality of fit significantly. 

The day-night asymmetry decreases with increase of Am^2 which is slightly disfavored 
by the data. Future more precise measurements of asymmetry will have stronger impact. 

2) Low energy observables, sensitive to pp- and Be- neutrino fluxes. They depend 
on the averaged vacuum oscillation probability. In particular, the Ge-production rate is 
proportional to 

Qoe oc cos^ ^13(1 - 0.5 sin^ 29). (16) 
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Notice that practically there is no dependence of these observables on Am^2- Since in the 
best fit point for zero 1-3 mixing Qtheor ~ Qexp, a significant increase of sin^ ^13 will lead to 
worsening of the fit (though small non-zero 1-3 mixing, sin^ ^^13 < 0.04, could be welcomed). 

3) KamLAND. The present best fit point is near the maximum of oscillation survival 
probability averaged over the atomic reactors. So, the increase of 1-3 mixing could be 
compensated by the decrease of 1-2 mixing. 

These features allow one to understand results of the data fit. 

We have performed analysis of the solar neutrino data for fixed value sin^ 6*13 = 0.067 
(fig. 4). The number of degrees of freedom is the same as in the 2u fit and we follow the 
procedure described in sect. 3 with survival probabilities given in (13). The best fit point 

Am^ = 11.0 • lO^^eV^ tan^ 9 = 0.38, /s = 1.03 (17) 

corresponds to /d.o.f. = 66.8/81. It is shifted to larger Am^ to satisfy condition (15). The 
shift to larger 1-2 mixing is disfavored by the Gallium experimental results (16). The allowed 
region is also shifted to larger Am^ and its size increases (also in the mixing direction). These 

results agree with the analysis in [15]. 

We also performed the combined analysis of the solar neutrino data and the KamLAND 
spectrum for sin^ ^13 = 0.067. In the best fit point we find xlnin/d-o.f. = 74.3/94 and 

Am^ = 7.1 • 10-W, tan^ 9 = 0.42, Jb = 1.03. (18) 

So, introduction of the 1-3 mixing slightly worsen the fit: Ax^ fa 1.9. The best fit point 
shifts to larger tan^^ to satisfy the condition (15). Significant shift to larger Am^ is not 

allowed by the KamLAND spectral data. At the same time, the 1-3 mixing improves a fit 
in the h-region. This region is accepted now at 90 % C.L. with respect to the best fit point 
(18). 

Notice however, that such an improvement is realized for values of sin^ ^^13 which are 
accepted by the CHOOZ and atmospheric data at the 3a level. So, inclusion of the CHOOZ 
result in the global fit changes a situation and no significant improvement of the fit in the 
h-LMA region occurs [18]. In fact, at sin^^i3 = 0.067 this region disappears even at the 3(7 
level. 

5 Physics of the solar neutrino conversion. 
5.1 Dynamics of conversion 

With new SNO results the h-LMA region is excluded or strongly disfavored and also sig- 
nificant deviation of 1-2 mixing from maximum is established. This essentially determines 
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both qualitatively and now quantitatively the physical picture of solar neutrino conversion 
(see also [26]). Different analyzes of data converge to 



Am^ = (6 - 8) • 10-^eV^ sin^ 6 = 0.28 - 0.30. (19) 

The non-zero 1-3 mixing may produce some small shift of the parameters. Since the physical 
picture is nearly the same for all points in these intervals, for further estimations we will 
take the values (8) as the reference point. 

We take the difference of the matter potentials for i/g and Va (active) according to the 
Standard Model: 

V = ^GppYjm^. (20) 

Here p is the matter density, is the number of electrons per nucleon, and is the 
nucleon mass. In fact, the latest experimental data allow to check the presence of such a 
potential in the model independent way. The extracted value of the potential is in agreement 
with (20) (within Icr) and K = is rejected at ~ 5.6(7 level [14]. 

For parameters (19) the neutrino evolution inside the Sun occurs in the highly adiabatic 
regime. It is convenient to write the averaged adiabatic survival probability as [3, 4] 

P = sin^ Q + cos cos^ (21) 

where is the mixing angle in the production point. Provided that there is no coherence 
lost (see discussion at the end of the section) , the depth of oscillations at the surface of the 
Sun equals [4] 

Ap = sin2^sin2^^. (22) 

So that the probabihty (being the oscillatory function of distance) is inscribed in the fol- 
lowing oscillation strip 

P± ^sin2^sin2^||^. (23) 

The oscillation length in matter, is smaller than the resonance oscillation length, 
Iji = Ij^/ sin 29, where li, is the vacuum oscillation length. So, for typical energy 10 MeV we 
find Im ~ 3.5 • 10^ cm. Along the solar radius, it!©, about Rq/Ith ~ 2 • 10^ oscillation length 
are obtained, and consequently, the oscillations are strongly averaged out. 

Dynamics of the effect depends on sin^ 9 and 

''^ ^~ Vo ~ 2EVo~ 2V2EGp{pYe)o' ^ ^ 

which is the ratio of the "kinetic" energy, Am'^/2E, and the potential energy in the neutrino 
production point, Vq. For a given energy, rj determines the relative contributions of the 
vacuum oscillation and the matter effect. 
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5.2 Two limits 



Depending on rj there are two limits: 

1. Matter dominance: r] <^ 1 which corresponds to — > oo. In this hmit the neutrino 
flavor evolution has a character of the non-oscillatory {Ap = 0) conversion with the survival 
probability [4] : 

Pnon—osc — sin 9. (25) 

In this case neutrinos produced far above the resonance density propagate to zero (small) 
density adiabatically. The initial mixing is strongly suppressed {9^ ~ 7r/2) and the propa- 
gating neutrino state practically coincides with the heaviest eigenstate: ^{t) — i'2m- At the 
exit from the Sun: u{t) = V2 [4]- 

2. Vacuum dominance: rj ^ 1 which corresponds to E ^ 0. Matter effects arc small. 
The flavor evolution has a character of vacuum oscillations, so P converges to the averaged 
oscillation probability 

P^Pva^^l- 0.5 sin^ 29, (26) 

and Ap ^ sm'^29. 

For the reference value of mixing (8) we flnd 

Pnon-osc = 0.281, P„„e = 0.596. (27) 



The resonance value of rj, for which P = 1/2, equals 

VR = = 2.28. (28) 

cos 29 

It marks the transition region between the two extreme cases. 

The boron neutrinos are produced in the region where the density pYg decreases from 
100 to 70 g/cc. The layer with maximum of emissivity (r = O.O43i?0) has the effective 
density pY^. ~ 93.4 g/cc [2]. For this density the value (28) corresponds to the resonance 
energy Eji = 2.2 MeV. Neutrinos with rj < rj^ {E > 2.2 MeV) are produced above the 
resonance (density) and then propagating to the surface of the Sun cross the resonance 
layer. Neutrinos with rj > rjp {E < 2.2 MeV) are produced below the resonance (density) 
and never cross the resonance (or cross the resonance layer twice depending on direction of 
propagation) . 

We can also define the "median" energy at which the probability equals the average of 
the two asymptotics (25) and (26): Pmed = ^-^iPvac + Pnon-osc)- According to (21), this 
corresponds to cos^ 9m = 0.5 cos^ 9 and is realized for 



Vmed 



COS 2^ 



2sm^^ ^ ^ 



Vl + sin^ 9, 



(29) 
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We find that rjmed = 1-43 and Emed = 3.54 MeV. 

The ratio CC/NC measured by SNO determines the survival probability averaged over 
the energy range above the SNO threshold for the CC events ~ 5 MeV): 

For the reference value (8) we get (P) = 0.322. This value is shghtly larger than the 
experimental result (1). 

So, both the experimental and theoretical values of CC/NC are rather close to Pnon-osc 
which means that at high energies {E > 5 MeV) the evolution of neutrino state is nearly 
non-oscillatory conversion. The difference 

CC/NC , CC/NC , , 

/ 1 = — ^-1 = 0.13 (31) 



■Pnon-osc Sin 9 



is due to 



• the averaged oscillation effect inside the Sun. In fact, for (pYf.)o ^ 93.4 g/cc and ener- 
gies relevant for the SNO CC signal, E = 14, 10, 6 MeV, we obtain 77 = 0.36, 0.50, 0.84 
correspondingly. These values of rj axe not small, though being smaller than the 
resonance value. Therefore one may expect a significant deviation of P from the 
asymptotic value. 

• the earth regeneration effect. 

5.3 Survived probability: adiabatic conversion versus oscillations 

The survival probability can be written as 

P = sin' e + APreg + APosc, (32) 

where APreg is the regeneration correction and APgsc is the oscillation correction. Both 
these corrections are positive. 

The regeneration effect, APreg, can be expressed in terms of the Day- Night asymmetry, 
Adn, as: 

^Preg ~ -^ND ' Pnon-osc ~ 

sin' 9 ■ And- (33) 
For the best fit point the asymmetry equals And — 3.0%, and consequently, APreg — 0.008. 
Using formula for the adiabatic conversion (21) we find: 

APo,c = cos2^cos'^^. (34) 
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Notice that cos^ 9^ gives the probabihty to find the eigenstate in the adiabaticaUy 
propagating neutrino state: |(z^im|^(^))P = cos^ = const. In the non-osciUatory hmit we 
would have cos^ 6'^ ~ and vit) = i'2m- The presence of second eigenstate, firm in ^ifyi 
leads to the interference effect, and consequently, to oscillations. 

In the limit of small 7] we find 

cos^ ^ sin^ 2^ r7^ (35) 

and therefore the correction to probability can be written as 

APo.c= ^cos29sui^29 {ri^ + 0{rf)). (36) 

Notice that the correction is quadratic in 77, and furthermore, it contains small pre-factor 
cos26'/4 ~ 0.1. It is for this reason the correction is rather small in spite of large values 
of r]. However, convergence of the series is determined by rj itself, and so, the corrections 
to the first order /^Posc are not small. Although (36) allows to understand the size of the 
correction, in our estimations we use exact expression for cos^ 9^^. In the limit of small r] 
the depth of oscillations 

Ap ^ sin^ 29 rj (37) 

decreases linearly, that is, slower than correction to the average value. 

In fig. 5 the averaged survival probability is shown as a function of the neutrino energy 
for different production points (different initial densities). The shadowed strips show the 
depth which oscillations would have at the surface of the Sun, provided that there is no loss 
of coherence. The average P converges to Pnon-osc (dashed line) with increase of energy and 
pYe. The decrease of oscillation depth with rj is much slower than convergence of P to sin^ 9: 
AP oc ?7^, Ap oc T]. The depth of oscillations increases with decrease of E approaching the 
vacuum value. Notice that even for the highest energies of the spectrum the conversion is 
not completely non-oscillatory, though P ^ sin^ 9. 

Using 1] = 0.50 which corresponds to a typical energy of the spectrum measured by SNO, 
E ^ 10 MeV, we get APosc ~ 0.030 (approximate formula (36) leads to APosc ~ 0.022). 
This together with the regeneration effect reproduces well the observed difference of P and 
Pnon-osc- The depth of oscillations for this set of parameters is rather large: Ap — 0.45. 

Thus, at energies relevant for the SNO CC events, the survival probability is about 
12% larger than the non-oscillatory probability. Oscillations give dominating effect in this 
difference. The regeneration contributes about 4%. 

At low energies {E < 2 MeV) the Earth regeneration (oc (Am^/EY) can be neglected 
and the probability is given by the vacuum oscillation formula with small matter corrections. 
For ?7 ^ 1 we can write: 

P ~ Pvac - 7^ COS 29 sin^ 29. (38) 
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For the beryllium neutrinos the effective density in the production region {pYe)o = 87 g/cc, 
and correspondingly, 7] = 6.28. Inserting this value of r) in (38) we find AP = —0.028 which 
is smaller than 5% of P^ac- 

For the pp-neutrinos the effective density in the region of the highest production rate is 
{pYe)o ~ 68 g/cc. At E = 0.4 MeV this gives rj = 17.3 and correction AP = —0.01. 

The results for different values of Am^ can be immediately obtained from the fig. 5 by 
rescaling of the energy: for Am'^ the probability at the energy E' equals P{E') = P{E) 
where E ^ E' ■ {Am^/Am'^). If, e.g., Am'^ = 14 • 10"^ eV^, we find from fig. 5, that 
P = 0.46, 0.40, 0.35 for E = 6, 10, 14 MeV. Notice that even for so large Am^ the 
probability is substantially lower than the vacuum value and for high energies the matter 
effect still dominates. 



5.4 Coherence loss 

Let us consider effect of the coherence loss of a neutrino state on the conversion picture. 
The loss of coherence suppresses depth of oscillations, so that the probability converges to 
the average value. 

The coherence length, Lcoh, can be estimated from the condition: 

/Lcoh 
dxAvm = S^, (39) 

where Avm = V2m — vim is the difference of the group velocities of the eigenstates, 5i, is the 
length of the wave packet and x is the distance. The group velocities equal: 

^im=^^: ^=1,2, (40) 

where Hi are the eigenvalues of the effective Hamiltonian of the neutrinos in matter. For 
the difference of the velocities we find [27]: 

Am^ r?3. — cos 29 , , , 

2p' ^772 -77^ cos 2^ + 1 

where = Am^ /2pV{x). 

Let us analyze the expression in (41). In vacuum, r]x ^ 00 and 

Am^ 

Avm = Av^ (42) 

With increase of density the difference Av^ decreases. In the resonance we find 

A 777 2 

A^;^ = ^sin2^. (43) 
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A^;^ changes the sign at r]x = cos 2^ which corresponds to density po = Pr/ cos^2d larger 
than the resonance density. For p > po the difference of velocities changes sign. At very 
high densities, rjx ^ 0, we have 

At;^ = -^cos2^. (44) 

Notice that in the resonance channel always |Av^| < \Av\. So, the matter suppresses the 
divergence of packets. Furthermore, due to change of sign of Avm the overlap of packets can 
be recovered. If neutrino propagates from high densities (above the resonance) the packets 
first diverge, but then below po the overlap and therefore the coherence can be restored 
again. 

For the smallest value of rjv ~ 0.36 in our case {E — 14 MeV, pYe — 94 g/cc) we get 
— — O.OSA^;. So, the divergence effect can be estimated using the vacuum value 

The sizes of wave packets are different for different components of the solar neutrino 
spectrum. Let us consider them in order. 

The lifetime of isolated nuclei, tb = 10^^ s, is very long and the size of the wave 
packet is determined by the average time between two consequent collisions, tcoii' 

Su = ctcoll- (46) 
Furthermore, the collisions with large enough momentum transfer Ap: 

Ap>^ (47) 

should be considered. For smaller Ap the coUision does not break coherence (shorten the 
packet). Using Eq. (46) and (47) we find the condition for S,^: 

5, = ctcoiiiAp > 2ti/8,). (48) 
The collision time can be estimated as 

tcoii = .\ , (49) 
a[/\p)vn 

where v is the velocity of colliding particles (protons, electrons, the helium atoms), n ~ 
p/rriN is the number density of these particles and a"(Ap) is the cross-section. The domi- 
nating effect is the Coulomb scattering, so that the cross-section can be written as 

a(Ap > = = p^. (50) 
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where in the lower approximation A ~ 7r{Za/v)^ corresponds to the Rutherford formula, 
where Z = 5 is the charge of Boron nuclei, a = e^/47r and v is the velocity of a colliding 
particle. In our case Za/v > 1 and the perturbation theory is violated. For estimations we 
take A ~ 0.1 — 1. Inserting (50) into (49) and then (48) we get 



which leads to 



nvA 



(51) 



(5^ = (1 -2) X 10-^cm. (52) 

This result is in agreement with results obtained by other methods [28, 29, 30] . The coher- 
ence length equals: 

Here Rq = 7 ■ 10^° cm is the solar radius. The matter effect leads to suppression of the 
divergence, and consequently, to larger coherence length. According to (51) for the observ- 
able part of the boron neutrino spectrum the effect of the coherence loss is small. Strong 
divergence of packets in the conversion region {x ~ O.IRq) occurs for £■ < 1.5 MeV. 

For the '^Be neutrinos the length of wave packet is determined by the production 
time [30], that is by the capture of the electron. The latter is given by the time an electron 
wave packet of the size Se crosses nucleus [30] : 

-I- (^^) 

Using the thermal velocity: ~ ^JZkT/me and the thermal wave length for the electron 
packet: 5e ~ 27r/meVe, we get from (54) for the neutrino wave packet length: 

A 27r 27r 

Numerically this gives ~ 6 • 10~* cm [30]. 

The Be neutrinos undergo nearly vacuum oscillations with small matter effect, so to 
calculate the coherence length we can use the formula (45): 

L,,h{Be) ~ O.Oli^o- (56) 

That is, the coherence is lost already in the production region. Notice that the oscillation 
length is even smaller (3 • 10^ cm) and Lcoh/L ~ 200. 

For the pp-neutrinos the coherence length is given by (45). The size of the packets is 
determined by the time of interaction which, in turn, is given by the wave packet of the 
colliding proton, similarly to the case of ^Be neutrinos. Notice that expression (55) does 
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not depend on a mass of colliding particle, and therefore it can be immediately applied for 
the pp-neutrinos. As a result, the neutrino wave packet has the same size (5,^ ~ 6 • 10~* cm. 
The coherence length is shorter due to smaller energies: 

L..,(pp)^ 0.002i?e(^;^)'. (57) 

Again the loss of coherence occurs already in the production region. The number of oscil- 
lation periods before complete averaging: Lcoh/lv < 100. 



6 Spectrum distortion at SNO 

Distortion of the energy spectrum is the generic consequence of the LMA MSW solution. 
As we discussed in the previous section, with decrease of energy the survival probability 
increases due to increase of the oscillation contribution. 

The energy spectrum of electrons has been calculated according to 

S{Teff)=[ f ^^%^R{Tr,Teff)\P''fBME)+P'''''<l>hep{E)]dEdTr , (58) 

where Tg// is the measured electron kinetic energy, is the real electron kinetic energy, E is 
the neutrino energy, (t)B{hep) is the Boron (hep) neutrino flux, da{Tr, E)/dTr is the differential 
cross section taken from [31], and R{Tr,Teff) is the SNO resolution function, and P^^^ 
are the survival probabilities averaged over the corresponding production regions. 

In fig. 6 we show the results of calculations of spectra for different values of Am^. The 
distortion due to oscillations which dominates at low energies is partly compensated by the 
regeneration effect at high energies. Thus, for the day signal one expects stronger upturn. 

According to fig. 6, the upturn is about 8 - 15%. We show also the SNO experimental 
points from the salt phase. A dependence of the distortion on Am^ is rather weak in the 
allowed region. Notice that in the low energy part the spectrum has a tendency to turn 
down in contrast to the expected effect. At < 7.5 MeV the points are systematically 
below the predicted rate. One should notice, however, that the experimental points include 
the statistical error only and it is not excluded that some systematics explains the observed 
result. 

The same effect - an absence of the upturn of the spectrum - is observed in the phase 
-I of the SNO experiment [5] . The spectral data agree well with the undistorted spectrum. 
It would be interesting to combine the results of both phases to improve statistics (and 
probably reduce the number of bins). Being confirmed, the fact of absence of the upturn 
or even a turn down at low energies can be explained by the effect of additional sterile 
neutrino 
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7 Next step 



Let us present predictions for the forthcoming experiments. 

1). KamLAND. More precise measurements of the rate and spectrum distortion are 
expected. That can further diminish uncertainty in the determination of Am^. In fig. 7 we 
show the contours of constant suppression of the KamLAND rate Rkl above 2.6 MeV: 

_ 7V(Am^tan^g) 

= TVo ' ^ ' 

where and Nq are the numbers of events with and without oscillations. As follows from 
the figure strengthening of the lower bound on Rkl will cut the allowed region from the 
side of small Am^ (~ 5 • 10~^ eV^) as well as large Am? (~ 9 • 10~^ eV^) and large mixings: 
tan^ 6 ~ 0.45. In contrast, strengthening of the upper bound on Rkl will disfavor the region 
of small mixings: tan^ 6 ~ 0.3. 

The spectrum distortion can be characterized by a relative suppression of rates at the 
high and low energies. We choose E — 4.3 MeV as the border fine [19], so that the 
interval (2.6 - 4.3) MeV contains 4 energy bins. Introducing the rates Rkl{< 4.3 MeV) and 
Rkl{> 4.3 MeV) we define the shape parameter as 

l-i?^i(<4.3MeV)- ^ ^ 

k does not depend on the normalization of spectrum and on the mixing angle in the 2v 
context. It increases with the oscillation suppression of signal at high energies, k > 1 
{k < 1) means stronger suppression at high (low) energies. The present KamLAND data 
give [19] 

^exp ^ 0.84l°:^2^ 1(7. (61) 

In fig. 7 we show the contours of constant shape parameter. According to this figure at the 
la level 

k'^ = 1.051°:^^ (62) 

in a very good agreement with (61). In the 1-LMA (la) region we find k = 0.5 — 2.3. So 
that even mild increase of statistics will influence the allowed range of Am?. 

In the h-LMA region the allowed interval, /c = 0.9 — 1.4, is narrower. So, if forthcoming 
measurements favor k < 0.9 or k > 1.4, the h-LMA region will be further discriminated. 

2). Precise measurements of the CC/NC ratio at SNO. In fig. 8 we show the contours of 
constant CC/NC ratio with finer grid than before. We find predictions for the best fit point 
and the 3a interval: 

^ = 0.32in?, (3a). (63) 
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3) . The day-night asymmetry at SNO. In fig. 8 we show also the contours of constant 
for the energy threshold 5.5 MeV. The best fit point prediction and the Scr bound equal 

A^^° = 3.0±0.8% , (la), A^^°<6% (3a). (64) 

4) . Germanium production rate. In the best fit point we predict Qce = 71 SNU. We show 
in fig. 9 the lines of constant Ge production rate with finer (than before) grid. 

For the Argon production we have Q^r = 2.96 SNU in the global b.f. point. We show 
in fig. 9 also the lines of constant Argon production rate. 



8 Conclusions 

1. The SNO-II fluxes have only slightly shifted the best fit point toward smaller Am^ and 
tan^ 9. The most important improvements consist however of the stronger upper bounds on 
Am^ and tan^ 9. Those imply that the h-LMA region is strongly disfavored and 1-2 mixing 
deviates substantially from the maximal one: 

(sin^ 9 - 0.5) ~ sin^ 9. (65) 

The 1-3 mixing does not change these results once the CHOOZ data are included in the 
analysis. 

2. These improvements in measurements of the oscillation parameters lead to a situation 
when physics of the solar neutrino conversion is essentially (and quantitatively) determined. 
In the high energy part of spectrum the averaged survival probability is close to the non- 
oscillatory one. For E > 5 MeV the effective P is about 12% higher than sin^ 9. Oscillations 
give the dominant contribution to this difference: ~ 8 — 14%, the rest is due to the Earth 
regeneration effect. In spite of smallness of difference (P — sin^ 9) which is proportional 
to 77^, the depth of oscillations is relatively large: Ap ~ 0.45 being proportional to rj. At 
low energies {E < 1 MeV), vacuum oscillations are the dominant process with the matter 
corrections to P below 5%. 

3. After the SNO salt results the errors of determination of the oscillation parameters 
become smaller than the values of parameters: 

5{Am'^) < Am^ ^(tan^ 9) < tan^ 9. (66) 

This means that the solar neutrino studies enter a stage of precision measurements. 

For the forthcoming measurements we predict about 8 - 15% upturn of the energy 
spectrum at SNO. If further measurements confirm the absence of the upturn hinted by the 
present data, some physics "beyond LMA" should be invoked. 

The CC/NC ratio is expected to be ^ 0.32, the Day-Night asymmetry: AfJl° ~ 3% 
(Tg// > 5.5 MeV), the spectrum shape parameter: A; = 1.0 -|- 0.8/ — 0.5. 
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Figure 1: The allowed regions in tan^ 9 — Am^ plane, from a combined analysis of the solar 
neutrino data at Icr, 90%, 95%, 99% and 3cr C.L.. The boron neutrino flux is treated as 
free parameter. The best fit point is marked by star. 
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Figure 2: The allowed regions in tan^ 9 — Am^ plane, from a combined analysis of the solar 
neutrino data and the KamLAND spectrum at la, 90%, 95%, 99% and 3cr C.L.. The boron 
neutrino flux is treated as free parameter. The best fit point is marked by star. 
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Figure 4: Three neutrino analysis with sin^ 6*13 = 0.067. The allowed regions in tan^ 6* — Am^ 
from a combined fit of the solar neutrino data at the la, 90%, 95%, 99% and 3a C.L.. 
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Figure 5: The averaged survival probability as a function of the neutrino energy for different 
production initial densities pYg- The oscillation strips (shadowed) show the depth of oscilla- 



tions at the surface of the Sun. 
is shown by the dashed line. 
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Figure 6: The energy spectrum of the CC reaction events. Shown is the distribution of 
events in the kinetic energy for different values of Am^ and tan^ 6 = 0.39. We show also 
the SNO experimental points from the salt phase. 
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Figure 7: The contours of constant rate suppression, Rkl (dotted) and the spectrum shape 
parameter k (dashed) at KamLAND. We show also the allowed regions of the oscillation 
parameters from the combined fit of the solar neutrino data and the KamLAND spectrum. 
The best fit point is indicated by a star. 
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Figure 8: Predictions for the CC/NC ratio and the Day-Night asymmetry at SNO. The 
dashed hnes are the hnes of constant CC/NC ratio (numbers at the curves) and the dotted 
hues show the hnes of constant (numbers at the curves in %). We show also the 

la and 3o" allowed regions of the oscillation parameters from the combined fit of the solar 
neutrino data and the KamLAND spectrum. The best fit point is indicated by a star. 
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Figure 9: Predictions for the Germanium and Argon production rates. The dotted hues are 
the hues of constant Germanium production rate, Qce, and the dashed hues show the hues 
of constant Argon production rate ^Qat (numbers at the curves in SNU), in the Am^ — tan^ Q 
plane. We show also the allowed regions of the oscillation parameters from the combined fit 
of the solar neutrino data and the KamLAND spectrum. The best fit point is indicated by 
a star. 
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